Non-equilibrium and proximity effects in superconductor-normal metal junctions 



CO 

o 

Oh 
< 






a 



T3 

O 

o 



> 

00 
00 
<N 

O 
CO 



1 



V. J. KauppiV'QH. Q. Nguyen, 1 and T. T. Heikkila 1 

1 0. V. Lounasmaa Laboratory, Aalto University, P.O. Box 15100, FI-00076 AALTO, Finland 

(Dated: April 5, 2013) 

We study the consequences of non-equilibrium heating and inverse proximity effect in normal metal 
- insulator - superconductor - insulator - normal metal (NISIN) junctions with a simple quasi-one- 
dimensional model. We especially focus on observables and parameter regions that are of interest 
in the design of SINIS coolers with quasiparticle traps. We present numerical results calculated 
by solving the Usadel equation and also present analytical approximations in two limiting cases: a 
short junction with a non-negligible resistance in both ends and a long junction with a transparent 
contact at one end. 



INTRODUCTION 

Electronic on-chip coolers offer a promising method for 
reaching below-lOOmK temperatures in variety of appli- 
cations including final stage of cooling in low temperature 
physics experiments as well as extremely sensitive radi- 
ation detectors used in space applications. A promising 
way to realize this type of cooler is by using a supercon- 
ductor (S) weakly coupled to the normal metal (N) which 
is to be cooled pQ. 

The operational principle of these NIS coolers is based 
on the presence of the energy gap in the superconductor 
density of states. When the bias voltage applied across 
the junction is adjusted correctly, only the hot quasipar- 
ticles can tunnel to the superconductor, thus cooling the 
normal metal. A more comprehensive review of the the- 
ory of NIS junctions can be found from Refs. [TJ [2] - 

The history of (SI) NIS coolers dates back to the 
nineties [2]. While the understanding of these systems 
has progressed, it has been understood that the main lim- 
itation to the cooling is often due to the non-equilibrium 
heating of the superconductor [3] . A popular solution to 
bypass this limitation is to use another normal metal, a 
"quasiparticle trap" , in contact with the superconductor 
to allow thermalization of the hot non-equilibrium quasi- 
particles in the superconductor [3l [4] . Later the same ef- 
fect was achieved by using magnetic fields [5], where the 
mechanism is essentially the same with normal metallic 
vortex cores acting as quasiparticle traps. Recently it has 
also been demonstrated that making the superconduc- 
tor wide close to the contact reduces the non-equilibrium 
heating [6]. However, this approach is not applicable in 
systems with wide junctions aiming at large cooling pow- 
ers. 

So far the non-ideal characteristics of NIS coolers has 
been analyzed with simplified thermal models, which as- 
sume the presence of a quasiequilibrium distribution in- 
side the superconductor, and do not include for example 
the inverse proximity effect. In this paper, our aim is 
to provide a microscopic description of the role of non- 
equilibrium effects in NIS coolers in the presence of a 
quasiparticle trap. We base our description on a quasi- 



one dimensional model and take into account the inverse 
proximity effect which also affects the cooling process in 
setups with well-coupled traps. Effects we do not take 
into account since they have been discussed elsewhere 
include the effect of the environment on the density of 
states [7 and coherent [8] and incoherent [9 Andreev 
effects. 



THE MODEL 

We model the NIS cooler with a quasi-one dimensional 
model shown in Fig. [I] The model consists of a normal 
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FIG. 1. Schematic picture of the quasi-one dimensional 
model. fFD,L/R are the Fermi-Dirac distribution functions 
in the left /right normal metals. 

metal island which is to be cooled (on the left) in equi- 
librium with some temperature T^, the normal metallic 
quasiparticle trap (on the right) in equilibrium with some 
temperature Tr and the superconducting layer of length 
L between the two. The superconducting layer is coupled 
to the two normal metals by insulating barriers with re- 
sistances Rl and Rr, similar to the geometry in [10 . 
The system is then biased with some voltage V. This 
model can also be used to study the NIS junction in the 
absence of the quasiparticle trap. In that case the length 
of the superconductor, L, must be set equal to the energy 
relaxation length of the superconductor. 

The effects we take into account in our calculations are 
1) the inverse proximity effect from the normal metals to 



the superconductor, 2) the non-equilibrium heating of the 
superconducting wire and 3) the electron-phonon inter- 
action in the normal metal island, which also adds two 
parameters to our theory: the electron-phonon coupling 
strength E multiplied by the volume E of the island and 
the phonon temperature T p h. 

To make the calculations tractable, we consider two 
type of coolers separately. The first type is one with a 
sufficiently long wire (in practice L > 3 — 4£, where £ is 
the coherence length of the superconductor) and a good 
contact between the quasiparticle trap and the supercon- 
ductor, i.e. Rr = 0. The second case considers a short 
superconducting wire (L < £) and arbitrary resistances 
on the two interfaces. 

Below, we switch to units where lengths are in the 
units of the superconductor coherence length, £, energies 
are in the units of the superconducting gap Ao at zero 
temperature and the resistances are in the units of the 
resistance of one coherence length long superconductor 
in the normal state, R^. We also set ks = h = e = 1 
so that the temperatures and voltages are also expressed 
in the units of energy. We retain SI units in the most 
important results. 

The starting point of our analysis is the Usadel equa- 
tion [11 which gives the quasiclassical Keldysh Greens 
functions of a superconductor in the dirty limit. The 
matrix Usadel equation can be further divided to two sep- 
arate equations [12] : the kinetic part which gives the dis- 
tribution functions in the superconductor and the spec- 
tral part which gives the spectral coefficients needed to 
solve the first one. The kinetic Usadel equation for the 
distribution functions of the superconductor reads 



d x (T> L (x)d x f L (x))=0 

d x (V T d x f T (x)) = f T (x)2TZ, 



(1) 



where we have divided the full distribution function of 
the superconductor to its symmetric and antisymmetric 
(or transversal and longitudal) parts in energy, f T and 
fL, respectively. The full distribution function can be 
obtained from these by / = \ (1 — Jl — /r). The co- 
ordinate and energy dependent spectral coefficients ap- 
pearing in Eq. u\l are the energy dependent spectral en- 
ergy and charge diffusion coefficients, T>l = cos 2 (Im#), 
T>t = cosh 2 (Re 0), and charge recombination coefficient, 
1Z = — Imsinh(#). These can be parametrized in terms 
of the pairing angle, 0(e, x), which can be solved from the 
spectral Usadel equation 



(Pi 
dx 2 



2iAcosh((9) - 2zesinh(0). 



(2) 



Generally, Eq. ([2| should also contain, and be accom- 
panied with another equation containing phase differ- 
ence, but for simplicity, we assume that the super current 
flowing through the superconductor is negligible so that 
phase gradient can be discarded [13]. The presence of 



the phase gradient would also add extra terms to the ki- 
netic equation. The three Usadel equations, two for the 
distribution functions and one for the pairing angle, are 
accompanied by six boundary conditions, which read [14] 

RL-Di(0)(d x fi(0)) = [fi(0) - f?{e,» L ,T L )] N s (0) 
RRDi(L)(d x fi(L)) = [f°(e, m ,T R ) - f^L)] N S (L), 
Rjd x 9\ x=0 ,L = sinh(6»(:r = 0,L)), 

(3) 
where Ns is the reduced density of states of the super- 
conductor and i G {L,T}. The reservoir distribution 
functions, ff are given by the symmetric and antisym- 
metric parts of the Fermi-Dirac distributions. It should 
be noted that due to the inverse proximity effect, the re- 
duced density of states in general differs from the bulk 
value calculated from the BCS theory, i.e. 



Nso = Re 



V^ 2 " 



A 2 



(4) 



After solving the three equations (HI) and pi) either 
numerically or analytically after some approximations, 
we can calculate the energy and charge current densities 
from the normal metal island by 



/CO 
deej L (x = 0) 
-oc 

/CO 
dej T (x = 0), 
-co 



(5) 



where ji = T>id x fi is the spectral current. The heat 
current from the island is then given by 



JQ=JE~ VLjc 



(6) 



In what follows, we disregard the self-consistency of 
the superconducting gap. The temperature dependence 
of the gap can be safely neglected since we are considering 
only temperatures Tl,Tr <C A. Also it has been shown 
[T5] that when Rr <C Rl, the non-equilibrium effects 
are not expected to change the gap considerably. This is 
also the limit at which the cooling is expected to work 
best and we may safely limit the calculations below to 
this limit. Furthermore, due to the proximity effect the 
gap close to the quasiparticle trap would be suppressed, 
which amounts to only renormalizing the length of the su- 
perconductor, leaving the qualitative results unchanged. 



DENSITY OF STATES 

As stated above, the most important effect from the 
proximity of the normal metal is to induce subgap states 
in the superconductor. To this end, we give an analytic 
expression for the proximity modified density of states in 
the two cases described above. 



Long wire, direct trap cooler 



Short wire, arbitrary resistances 



In this case, we have two contributions to the density 
of states that can be separated. First is the modification 
due to the right (good) interface and the second is the 
modification due to the left (high resistance) interface. 
To calculate the former one, we solve the spectral Usadel 
equation (J2| in the limit of a long wire and Rl ^> 1, 
Rr = 0. The result is 



exp (-(L - x)V2aj tanh(<9 5 /4) 



0{x) ~ 0s — 4 arctanh 

(7) 

where 0s = arctanh (A/e) is the pairing angle of a bulk 
superconductor and a = VA 2 — e 2 . This solution de- 
scribes a semi-infinite wire with a good contact to a nor- 
mal metal reservoir at x = L. The density of states can 
now be calculated from N$ = Re cosh 0. With the pairing 
angle (|7|) this becomes 



N s = N so + Re 



e(l + 6/3 2 +/3 4 )-4A(/3 + /3 3 ) 



V^ 



A 2 (l-/3 2 ) 



2\2 



(8) 



where (3 = exp(— L\r2a) tanh(0£/4). By expanding the 
second term for small /?, we get outside the region |A — 

e\ < 1/(2AL 4 ) 



When L <C 1, we can calculate the pairing angle as 
follows. We approximate the derivative in Eq. Q with 
a difference, d 2 x « (d x 0(x = L) - d x 0(x = 0))/L. Then 
we use the boundary conditions and assume that does 
not depend on x. The solution becomes 
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FIG. 2. (Color online) Density of states at the left end of the 
wire as a function of energy with Rl = 10000 and Rr = 0. 



N s = N Sfi - Re 



4Ae 



-LV2a 



tanh(6» s /4) 



Ve 



A 2 



(9) 



At zero energy, this becomes Ns(e = 0) « 1.7e L ^ 2 . 

The latter contribution, i.e., the contribution due to 
the finite (but large) resistance in the left interface can 
be calculated by expanding around its bulk value and 
then solving the spectral Usadel equation. The density 
of states from this becomes 



where R = RlRr/(Rl + Rr) is the parallel resistance 
of the two barriers and 7 = A^R^/(2LR) is the Dynes 
factor. The density of states from this becomes of the 
Dynes form 



N s = Re 



(e + i 7 ) 2 



V(e + z 7 ) 2 -A 2 



(12) 



It is worth noting that in this limit, the Dynes parameter 
is independent of energy. 



N s (x = 0) w Nso + Re 



v / 2i^L(A 2 -e 2 ) 5 / 4 ' 



(10) 



When the deviation of the pairing angle from the bulk 
value is small, these two contributions can be added to 
give the total density of states in this geometry. The 
rounding of the density of states for different lengths of 
the wire using an exact numerical solution for some pa- 
rameter values is shown in Fig. [2j 



HEAT CURRENT 

The heat current from the normal metal in an ideal 
NIS junction (i.e., without taking into account the non- 
equilibrium heating or the proximity effect), assuming 
the superconductor and the normal metal are both in 
equilibrium with some temperatures Tl, Tr and distribu- 
tion functions given by Fermi-Dirac distributions, JfDi 
reads 16 



Jq IS = jrf d <t ~ V ) N S0 (f FD (e,V,T L ) - f FD (e,0,T R )) 



R L V 2A3 



,53/2 



V 



T L 



V 



T L 



-9l/2 



V 



T L 



A 2 
Rl 



2T r tt c _ a/Tr 



(13) 



where the approximation holds when V + Tl < A. We 
now proceed to considering the two additional modifi- 
cations to this - non-equilibrium and proximity effect - 
separately as corrections to this formula in the two ge- 
ometries descibed above. 



Inverse proximity effect 

For a moment we forget the non-equilibrium effects 
and write the heat current as in the ideal NIS case, but 
with the proximity modified density of states, i.e. 



jQ E =\Jde(e-^ L )N s 



(14) 



x Ufd{£,V>l,Tl) ~ fFD(e,fJ>R,T R )) , 



where R = Rl + Rr and \±l — jur = V. For voltages 
V+T^ <C A, we can expand the density of states at e « 
and calculate the integral. Formally the heat current 
density can be written as 



ape _ -NIS 
Jq - Jq 



■s$ E , 



(15) 



where the first term is the ideal NIS current given above 
and the second term is a correction to this. 

In the first type of a geometry with a direct trap and a 
long wire using the zero-energy expansion the correction 
becomes 



H 



PE 



2A 2 
''^Rl 



tan 



3A 2 






-f\/2 



Tl) 



eV 



(16) 



In the second type of a geometry, i.e. L <C 1, we 
get after identical considerations, in the relevant limit 
Rl ^> Rr 



Sj 



PE 



A7 

~^R 



6A 2 



{Tl 



Tl 



eVY 



(17) 



Next we proceed to consider the effect of the non- 
equilibrium heating. 



we calculate them in the first type of a geometry with 
a long wire and Rr = 0. We approximate the spectral 
coefficients by their bulk values and solve Eqs. ([l]). For 
the longitudal spectral current density at the interface of 
the island and the superconductor we get 



RlJi 



NsoSf L (e) 



1 



L 
Rl 



V l N s 



and for the transversal current density 



RlJt = 



NsoMe,V) 



N s tanh(^1Z/V T L) ' 
T> T R L 



(18) 



(19) 



1 



Analysis of these shows that the modification to the 
transversal current density vanishes when e = A and 
is proportional to ~ L/Rl otherwise and can thus be 
neglected since L/Rl < 1 for a typical setup. The cor- 



rection to Jl in expression (18), however, while also pro- 
portional to ~ L/Rl, diverges when e « A and thus 
needs to be taken into account. To calculate the correc- 
tion, we expand the energy mode of the heat current, i.e. 



integral of e times Eq. (18), in the limit of small L/Rl- 
The result becomes in the first order 



Rl Je « RlJe 



NIS 



Rl63 N e EQ , 



(20) 



where 



R L 5j N E EQ 



LRf: 

'Rti 



/>oo 
JA+u 



A 2 



A/ L (e). (21) 



Here we have also introduced a cutoff energy, ^, to cut 
off the logarithmic divergence of the integral at e — >• A. 
In the limit V + 2T L /r < A we can approximate the 
distribution functions with exponential functions and get 
as a result 



5j N E EQ 



A 2 LR^ 
^2Rtt ( 



3 T L log 



k B T L 



(22) 



where 7 « 1.78 is the Euler's constant. The cutoff is 
chosen from the condition (L/Rl)Ns ~ 1, which gives 
v~/\L 2 Rl/(2R 2 L i 2 ). 

As for the second type of geometry, i.e., L«l and 
Rr > 0> the same expression (22) holds, provided that 
Rl > Rr. 



Non-equilibrium heating 



Total heat current 



Now we need the non-equilibrium expressions for the 
heat current that can be solved from Eq. (fTl). First 



In the limit where the two corrections to the ideal NIS 
current are both small, we can get the total current from 



the NIS current by adding the two corrections together, 
i.e. 



JQ 



j% IS 



■PE 



sj^ Q 



(23) 



We first consider the total heat current for the setup 
with a long wire and a good contact between the trap and 
the superconductor. Since the proximity effect inhibits 
the cooling in short superconductors and non-equilibrium 
effect inhibits it in long superconductors, there exists a 
local maximum of the cooling as a function of the length 
of the superconductor. This behavior is shown in Fig. [3] 
Along with this analytic approximation we have plotted 
the heat current calculated by numerically solving the 
Usadel equation. The main deviations from the numer- 
ical data are due to the breakdown of the semi-infinite 
approximation at small lengths of the superconductor. 
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FIG. 3. (Color online) Heat current as a function of the length 
of the superconductor with eV — 0.7Ao and UbTl — ksTR — 
O.lAo- The dashed line is the exact numerical result and 
the solid line is approximation (23). The black solid line is 
the ideal NIS current. The upper lines closer to the ideal 
NIS current (red) are for Rl/R^ — 10000 and the lower lines 
(blue) are for R L /R^ = 1000. 

The optimal length of the wire is found from expression 



(23). It is 



L, 



opt 



V2 



log 



4v / 2tan(|)(^(rE-rl)-F 2 )i? L 



l°g (if) ^ 



(24) 

This expression is expected to work when the assump- 
tions made in the calculation of the heat current correc- 
tions above hold. The optimal length and the heat cur- 
rent at that length as calculated from the full numerics 
are shown in Fig. [4j For large Rl, the behavior corre- 
sponds roughly to Eq. ([24|, albeit there is a difference of 
L ~ 1 between the two, which results from our approxi- 
mation to use the bulk values for the spectral quantities 




FIG. 4. (Color online) Numerically calculated optimal length 
(upper figure) and the heat current at that length (lower 
figure) as a function of the interface resistance at eV = 
A - 0.66k B T and T L = T R = T. In the upper figure, we 
have also plotted the analytic approximation (24) (dashed 
line) shifted by a constant to account for the discrepancy dis- 
cussed in the text. 



Temperature of the island 

Next we consider the minimum temperature the nor- 
mal metal can be cooled into. This temperature mini- 
mum can be found by solving for T^ the heat balance 
equation 



JQ = 0. 



(25) 



For an ideal NIS junction, the limit is found to be 
T min _ i.42A(27rT i? /A) 1/3 exp(-2A/3T i? ). In practice, 
however, the non-idealities limit the minimum tempera- 
ture. We now consider these limitations. 

First we consider the long- wire setup. If the wire is 
long enough, the proximity effect does not matter and 
we can take into the heat current only terms Jq IS and 

SJq . We also assume that the electron-phonon cou- 
pling is negligible and instead analyze it separately below. 
The minimum temperature becomes 



j A 
k B 



A 



k B T R 



LRz 



2 t2 



B log 



k B T L Rji 
AL 2 RJ 



-i 2/3 



(26) 



where A and B are approximately constants at the opti- 
mal bias voltage when the ratio (A — eV)/kBTL is fixed. 
A comes from the NIS current expression and is approx- 
imately A « 4.25 at the optimal bias voltage. B is given 
by the non-equilibrium correction to the current and its 
numerically found value is B w 1.7 at the optimal bias 



voltage. From expression (26) we can see that at very 



low environment temperatures Tr, the minimum temper- 
ature is limited by ~ {L\og{RL / L) / Rl) • Numerical 
simulations showing this behavior are shown in Fig. [5] 
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FIG. 5. (Color online) Minimum temperature of the normal 
metal island as a function of the temperature of the trap with 



L = 13f and eV « eV c 



opt 



A - 0.66k B T L . For large T R , it 



is limited by the first term in Eq. ( 26 ) whereas for low Tr it 



saturates to the value given by the second term. 

In the case of a very large resistance Rl, the limita- 
tions to the temperature are no longer due to the non- 
equilibrium heating but due to the electron-phonon cou- 
pling and thus this effect must be added to the heat cur- 
rent. The electron- acoustic phonon heat current is given 
by Q] 



JQ 



-ph 



a(T, 



T p h) 



(27) 



where a = Ef2e 2 i2^A 3 /k b B is the dimensionless electron- 
phonon coupling constant. In this case, as Tr —> 0, the 
minimum temperature is given by 



Tmiin 
L 



A 
k B 



2R 2 L z 2 n 2 



nA^k™ 



1/3 



kBTph\ 
A J 



10/3 



(28) 



In the short junction setup we have different kind of a 
competition between the non-equilibrium and the inverse 
proximity effect. In the relevant limit Rl ^> Rr the 
proximity effect always dominates the non-equilibrium 
term in the heat current. This can be seen by comparing 
the prefactors of the two and noticing that 



R5j P Q E 



OC 7 : 



1 ^ L D r -NEQ 



(29) 



However, if we make the right resistance larger, at some 
point the non-equilibrium heating starts to matter and 
then it limits the cooling. Thus there must be an opti- 
mum value of Rl/Rr which maximizes the cooling. We 
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FIG. 6. (Color online) Equilibrium temperature of the normal 
island as function of Rl/Rr- Blue, red and green lines are 
L/£ = 0.5, 1, 2 respectively. Solid lines are for R L /R^ = 10000 
and dashed lines are for Rl/R$ = 1000. The temperature of 
the trap is k B T R /A = 0.1. 



show this behaviour in Fig. [6j These considerations are 
valid until we increase the total resistance of the junc- 
tion enough so that the electron-phonon coupling starts 
to dominate. In this case, the minimum temperature is 



again given by Eq. (28). 



CHARGE CURRENT 

We repeat similar analysis as above for charge current. 
For an ideal NIS junction the charge current is given by 

RI NIS = J deN S0 (f FD (e,V,T L ) - f FD (e,0,T R )). (30) 

For bias voltages V <C A this becomes 



NIS .v<* V2n Ak B T L A/kBTL / eV^\ 

e \k B T L ) 

(31) 
For bias voltages V ^> A the NIS charge current reduces 
to the Ohms law plus a correction that vanishes at high 
voltages: 

1 A 2 
V -^^. (32) 



NIS 



RI - ' 2e*V 

Using similar methods as above to calculate the heat 
current, we can calculate the corrections to NIS current 
due to non-equilibrium effects and the inverse proximity 
effects. We first calculate the long junction case. The 
non-equilibrium correction becomes (to the first order in 
L/R L ) 



( k B T L -(A-eV)/k B T L LRj P V « A 



eV» A 



(33) 



For V ^> A the correction due to the proximity effect 
vanishes, but for V <C A it becomes 

R L 5I PE « 4 tan (!) e~ L ^^V. (34) 

For the short junction setup, we calculate only the 
proximity correction. For V <C A it is given by 



R L 6I PE *^[def(e-» L ) = 



VZRJ 
LR L R R ' 



(35) 



We note that unlike the heat current (22), the non- 



equilibrium correction to the charge current at low bias 
voltages is proportional to the temperature, making it 
less important than in the case of the heat current. This 
can be traced back to the fact that the charge current is 
given by only the transversal part of the spectral current 
in which, as we noted above, the coefficient in front of 
the small parameter L/Rl does not diverge at e = A. 



CONCLUSIONS 

We have considered the heat and charge transport in 
NISIN junctions. We derived analytical estimates for the 
currents and temperatures of the normal metal island in 
two types of setups relevant for NIS cooler experiments: 
short junction with non-transparent contacts and long 
junction with a transparent contact between the super- 
conductor and the quasiparticle trap. We discussed three 
different effects on the heat current: the inverse prox- 
imity effect to the superconductor, the non-equilibrium 
heating of the superconductor and the electron-phonon 
coupling. 
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